J Stat Phys manuscript No. 

(will be inserted by the editor) 



Ladislav Samaj 

Charge and Current Sum Rules in 
Quantum Media Coupled to Radiation 



Received: / Accepted: 



Abstract This paper concerns the equilibrium bulk charge and current den- 
sity correlation functions in quantum media, conductors and dielectrics, fully 
coupled to the radiation (the retarded regime). A sequence of static and time- 
dependent sum rules, which fix the values of certain moments of the charge 
and current density correlation functions, is obtained by using Rytov's fluc- 
tuational electrodynamics. A technique is developed to extract the classical 
and purely quantum-mechanical parts of these sum rules. The sum rules are 
critically tested in the classical limit and on the jellium model. A compari- 
son is made with microscopic approaches to systems of particles interacting 
through Coulomb forces only (the non-retarded regime). In contrast with 
microscopic results, the current-current density correlation function is found 
to be integrable in space, in both classical and quantum regimes. 
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1 Introduction 

This paper is about the equilibrium statistical mechanics of infinite (bulk) 
media, conductors and dielectrics, fully coupled to the radiated electromag- 
netic (EM) field. The crucial problem of statistical analysis of such models is 
to determine how fluctuations of various statistical quantities, like charge and 
current densities, scalar and vector potentials, electric and magnetic fields, 
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etc., around their mean values are correlated in space and time. A special at- 
tention is devoted to the behavior of correlation functions at asymptotically 
large distances and to the so-called sum rules, which fix the values of certain 
moments of the correlation functions. 

The studied models are composed of spinless charged particles which are 
quantum, but non-relativistic (i.e., they obey the Schrodinger equation and 
not the Dirac particle- antiparticle formalism); according to the correspon- 
dence principle, a quantum system admits the classical description in the 
high-temperature limit. On the other hand, the interaction of charged parti- 
cles via the radiated EM field can be considered either non-relativistic (non- 
retarded) or relativistic (retarded). In the non- retarded regime, the speed of 
light c is taken to be infinitely large, c = oo, ignoring in this way magnetic 
forces acting on the particles; in other words, charges interact pairwisely only 
via instantaneous Coulomb potentials. In the retarded regime, c is assumed 
finite and the particles arc fully coupled to both electric and magnetic (trans- 
verse) parts of the radiated field. The ultimate aim is to describe quantum 
particles in the retarded regime, but due to serious technical difficulties one 
usually starts with the simpler non-retarded regime. 

Two complementary approaches exist in the theoretical study of charged 
systems. The microscopic one is based on the explicit solution of models de- 
fined by their microscopic Hamiltonians. The macroscopic one is based on the 
assumption of validity of macroscopic electrodynamics. Microscopic descrip- 
tion is complicated and laborious, but if some explicit results are available 
their value is considerable and, e.g., they can reveal a restricted applicability 
of macroscopic methods. 

Microscopic approaches are usually restricted to the non-retarded regime. 
A particular interest was devoted in the past to classical and quantum models 
of conductors. Classical particles interacting by the instantaneous Coulomb 
potential have been treated in the pioneering work by Debye and Hiickel [1] . 
An exponential shielding of the Coulomb potential, typical for the classical 
regime, was observed. One of the first classical analyses of a system of charged 
particles fully coupled to the radiation, based on the linearized Vlasov equa- 
tion and Maxwell's equations, was performed by Felderhof 2,3J, see also the 
monograph [4j. The existence of general sum rules for certain moments of 
the classical charge density correlation functions was established by Stillinger 
and Lovett [5]. The quantum sum rules, static and time-dependent, were the 
subject of numerous studies fOj^T^SuSslQl- A treatment of quantum Coulomb 
fiuids, based on the use of the Feynman-Kac path integral representation 
of the thermal Gibbs weight [TT], indicates that the charge-charge density 
correlation functions exhibit a long-range 1/r^ decay (see review [12] )■ A mi- 
croscopic analysis of conductor systems when transverse EM interactions are 
added to the Coulomb ones, was done recently p!3llT4] : a conceptual problem 
related to the use of classical radiation was pointed out in [TS]. Microscopic 
models of dielectrics are less developed; we mention a fluid model of quantized 
polarizable particles p!6lfT7| with the static dipole-dipole interaction. 

Macroscopic phenomenological approaches usually allow us to predict, 
with less effort, basic features of relatively complicated complex physical 
systems in the retarded regime. A macroscopic theory of equilibrium ther- 
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mal fluctuations of the EM field in conductors and dielectrics was proposed 
by Rytov [T8 1 IT9 1 I20] , see also the Landau and Lifshitz Course of Theoretical 
Physics [21] and the book [22]. For bulk systems, an interesting phenomenon 
occurs in two-point correlation functions of microscopic quantities like the 
induced electric potential or field: the presence of transverse interactions be- 
yond Coulomb cancel out some (spatially nonintegrable) long-ranged terms 
observed in the non-retarded regime [T5ll21j . Main applications of Rytov's 
theory were related to inhomogeneous situations. In the famous Lifshitz pa- 
per j23j, Rytov's theory was used in the context of the Casimir attraction 
between parallel dielectric plates. Another application appeared in the field 
of thermally excited surface EM waves [24]. Recently, the theory was used 
in the study of retardation effects on the long-range surface charge density 
correlation function between two distinct media [251126]. 

Macroscopic approaches, being essentially of mean-field type, are ex- 
pected to provide reliable results for the leading terms in the asymptotic 
long-wavelength behavior of correlations. Those results can be recast as sum 
rules for moments of correlations integrated in real space, or as asymptotic 
large-distance behaviors. As concerns these asymptotic behaviors, macro- 
scopic mean-field teories fail when fluctuations become crucial. The failure 
of mean-field predictions is well illustrated by the existence of 1 jr^ tails in 
static charge-charge correlations of quantum Coulomb systems, as explained 
in 2"^; see also the exact calculation of the tail for hydrogen plasma at low 
densities |28j . Another implicit assumption for the validity of macroscopic 
approaches is that correlations decay sufficiently fast. It is then legitimate to 
perform local averages for EM fields slowly varying in space, which reflects 
itself in the introduction of the frequency dependent dielectric function e(w) 
in Maxwell's equations. 

The present paper concerns sum rules for the bulk charge and current 
density correlation functions in quantum media, conductors and dielectrics, 
fully coupled to the radiation. Although this subject (in the non-retarded 
version) was in the center of interest of microscopic theories in the past, 
it remained practically unheeded by macroscopic approaches. Here, a se- 
quence of static or time-dependent sum rules, known or new, is obtained by 
Rytov's fluctuational electrodynamics. A technique is developed to extract 
the classical and purely quantum-mechanical parts of these sum rules. The 
sum rules are critically tested on a jellium model of conductors. A compari- 
son is made with microscopic approaches to systems of particles interacting 
through Coulomb forces only; in contrast to these microscopic approaches, 
the current-current density correlation function is found to be integrable in 
space, in both classical and quantum regimes. 

The paper is organized as follows. In Sect. 2, we review briefly Rytov's 
theory of EM field fluctuations and present basic expressions for the charge 
and current densities. The sum rules for the charge-charge, charge-current 
and current-current density correlation functions are derived and discussed 
in Sects. 3, 4 and 5, respectively. Sect. 6 deals with problematic points of both 
the Rytov fluctuational theory and microscopic approaches. Sect. 7 is devoted 
to an interesting comparison of the present results to those of Felderhof [21 
[5, in the classical limit. Sect. 8 is the Conclusion. 
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2 Macroscopic fluctuational formalism 

We consider the (3 + l)-dimensional space, defined by Euclidean vectors 
r = (xi, X2, X3) and time t. The two-point functions considered in this paper 
will be translationally invariant in both Euclidean space and time, so we can 
use the spectral (Fourier) representation 

fit,r;t',r')^f{t-t\r-r')= J^J ^e-M'-*')+''i(-'-')/(^, q), 

(2.1) 

where lo denotes the frequency and q = (91, 92, 93) the wave vector. The 
physical system of interest is a medium (conductor or dielectric) and an EM 
field present in it, which are in thermal equilibrium. 

The medium is composed of moving charged particles which are assumed 
to be non-relativistic. In the long-wavelength scale much larger than the 
mean interparticle distance, the macroscopic characteristic of the medium 
is the frequency dependent dielectric function e(w). We shall assume that 
the matter has no magnetic structure, i.e. it is not magnetoactive, and per- 
meability /i = 1. The medium is coupled to the EM field generated by the 
particles. The classical EM field is described by the scalar potential (j){t, r) 
and the vector potential A{t, r) with component Aj{t, r) (j — 1, 2, 3). In the 
considered Weyl gauge = 0, the microscopic electric and magnetic fields 
are given by 



c 



dt 



B = curl A, Bj = J2 ""jkij^Ai (2.3) 
k,i "^^ 

with c being the velocity of light and Cju being the unit antisymmetric pseu- 
dotensor. The elementary excitations of the quantized EM field are described 
by the photon operators of the vector potential components Aj which are 
self-conjugate Bose operators. The retarded photon Green function tensor D 
is defined by 

{A,(t,v)Ak{t' y) - Ak(t' y)A,{t,r)), t > t', 
iD,k{t-t';r,r') = { (2.4) 
0, t < t'. 

Here, Aj{t, r) denotes the vector-potential operator in the Heisenberg picture 
and the angular brackets represent the equilibrium averaging at temperature 
T. For non-magnetoactive media, the Green function tensor possesses the 
symmetry 

Djk{t - t';r,r') ^ Dk,{t - t';r' ,r). (2.5) 

The EM fields are random variables which fluctuate around their mean 
values governed by Maxwell's equations. To describe fluctuations of the EM 
fields, Rytov p^fT9 l [20 l [2T] studied the effect of a weak classical current j{t, r) 
due to the random particle motion in the medium. This current acts as an 
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"external force" on the vector-potential operator in the interaction Hamilto- 
nian 

Wint(t) = -^y"dVj(t,r).A. (2.6) 

The mean values of the components of the vector-potential operator can be 
thus expressed in terms of the Green function (|2.4p by using Kubo's linear 
response in current, 



~ J d\' j dt'Y,D,,{t-t';r,r')Mt',r'). (2.7) 



The mean value A satisfies macroscopic Maxwell's equations due to the clas- 
sical current j. This fact automatically implies a set of differential equations 
of dyadic type fulfilled by the Green function tensor. In particular, in the 
frequency Fourier space, we have 



E 



dxjdxi c 



Dikiu;;r,r') = -inhSjkSir - r'). (2.. 



This set of equations has to be supplemented by the obvious boundary con- 
dition of regularity I?;fc(w; r, r') for the distance |r — r'| going to infinity. 
In the case of a spatially homogeneous infinite medium, Djk{t — t';r,r') = 
Djk{t — t'; r — r'), under the Fourier transform with respect to the difference 
r — r' Eq. (|2.8p reduces to a set of algebraic equations, whose solution reads 

oj'^eyuj) oj^^eluj) — [cqy \ q'^ ' 



Here, q = |q|. The first/second term on the right-hand-side (r.h.s.) of this 
expression represents the longitudinal/transverse component of the Green 
function. Formula (|2.9p does reduce to its well-known form in the vacuum 
when one sets e{uj) — 1. 

Eq. (|2.7p defines —Djk{uj;r,r')/{h,c'^) as the tensor of susceptibilities cor- 
responding to the random variable A/c. The fluctuation-dissipation theorem 
tells us that the fluctuations of random variables are expressible in terms 
of the corresponding generalized susceptibilities. For the present symmetry 
(|2.5p . the theorem implies 



(A.AkZ^^ = -coth ( ^ ) lmDjk{to,q), (2.10) 



where f3 = l/{kBT) is the inverse temperature and {AjAk)lj q is the Fourier 
transform of the symmetrized correlation function 

{A,{t,v)Ak{t' ,v')Y ^ i(4(t,r)i,.(t',r')+ife(i',r')i,(t,r)). (2.11) 

We note that the symmetrization (|2.1ip is needless in the classical case. 

The fluctuation formula (|2.10p enables us to calculate the correlations of 
arbitrary statistical quantities provided that they are expressible in terms 
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of the components of the vector-field potential, within the classical EM field 
and in the gauge defined by Eqs. (|2.2p and p.3p . Let a scalar quantity u be 
expressible as u{t,Y) = U • A(<, r) = UjAj{t,T), where Uj are operators 
acting on t and r variables. Within the spectral representation with a single 
frequency and wave vector, f{t,r) — e~'"*+"i '"/(ti;, q), this relation takes an 
algebraic form u{LO,q) = Uj{uj,q)Aj{LO,q). It follows from the definition 
(|2.ip that the Fourier transform of the symmetrized correlation function of 
statistical quantities u and v, {u{t,r)v{t' ,r'))^ , is then determined by 

Mi,q = (u(t^,q)v(-t^,-q))' 

= J2 q)^fe (-^, -q) (Aj {to, q)Ak i-uJ, -q)) ^ , (2.12) 

where (A^ (w, q)Afc(-w, -q))^ = {AjAk)t,q given by (l2l0l) . 

In this paper, we shall concentrate on the volume charge density p and 
the electric current density j. The charge density is expressible as 



p(t,r).-LdivE = --i-^ 



d d 



An ' Attc ^ dx-j dt 



Aj{t,v) 



(2.13) 



or, in the spectral representation, as 

Anc 



p(w,q) = -^^^(IjAjiuj.q). 



(2.14) 



The electric current density is expressible in terms of the electric and mag- 
netic fields as follows 



j(t,r) = ^curlB(t,r)-^|E(t,r). 



(2.15) 



Substituting here the gauge relations (|2.2p and (|2.3p . and using the summa- 
tion formula for the elements of the unit antisymmetric pseudotensor 



(2.16) 



we obtain 



E 

L I 



92 



dxidxk 



or, in the spectral representation. 



1 92 



(2.17) 



jfc(w,q) = 



An 



-(lk^qiAi{oj,c{) 



q 



A-(t^,q) 



(2.18) 



It is a simple task to verify that both versions of the continuity equation 
d 



dt 



p(i,r) +divj(t,r) = 0, 



-upiLU,q) + J2lkjki^,q) ^0 (2.19) 
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are satisfied. 

Henceforth, in order to simplify the notation we set t' = and r' — 0, i.e. 
the time and position differences between the two points in the correlation 
functions will be t and r, respectively. 



3 Charge-charge density correlations 

We start with the charge-charge density correlation function. Using the spec- 
tral representation (|2.14[) it is easy to show that 

/^(PP>i. - ^ (£)'E'^.'^^-(^.^0^„ = -'Z^i^Im^, (3.1) 

where 

,(.).^cothr^V (3.2) 



Note that g{uj) > 1; the equality g{uj) ~ 1 takes place in the limit PTiu! 0. 
The fact that there is only one term of the order in the Fourier represen- 
tation p.ip is related to the applicability of the fluctuational theory in the 
long- wavelength limit q * (characteristic length) (characteristic length is 
usually the mean interparticle distance); an exact formula for P{pp)^ q, valid 
for any value of q, would contain also higher powers of q^. 

The inverse Fourier transform of p.ip leads to the zeroth-moment (neu- 
trality) condition 



d^r/3(p(i,r)p(0,O))« = (3.3) 
and the second-moment condition 



1 / d^rr'(3{pit,r)p{0,O)r = H dco exp(-ic.i)|Mlm ^ 
3 J t{u} 



(3.4) 



These sum rules do not depend on c and so they are the same in both non- 
retarded and retarded regimes. Note that the sum rule (|3.4p is not universal 
in the sense that its r.h.s. is a complicated function of the temperature and 
the medium characteristics. 



3.1 Static case 



In the static i = case, using the notation /5(0,r) = p{r), the general result 
can be rewritten in the form 



dVr"/3(p(r)/5(0))'^ 



— Im/(cj , / w = — - 
ui 27r^ 



1 



(3.5) 

note that the subtraction of 1 from l/e(u;) has no effect due to the operation 
Im. 



Im CO 




E C, Re CO 

Fig. 1 The contour in the complex frequency plane. 



The integral on the r.h.s. of Eq. (|3.5p can be formally evaluated by us- 
ing integration techniques and the general analytic properties of dielectric 
functions in the complex frequency upper half-plane, in close analogy with 
the treatment of the surface charge density correlation function in Sect. IV 
of Ref. PB]. We shall not repeat the derivation procedure and only write 
down the final result. Although our /-function defined in p.Sp differs from 
that considered in [26 , it possesses the necessary symmetry f*{w) = f{—uj) 
for real uj and goes to zero in the asymptotic limit — s- oo due to the 
asymptotic relation pMSU] 



(3.6) 



Here, tUp is the plasma frequency of the medium composed of species a (elec- 
trons and ions) with the number density n^, charge Ca and mass m^; since 
the masses of ions is much bigger than the mass of electron, it is usually suffi- 
cient to approximate the plasma frequency of the medium by aj^ = Anne'^ /m 
with n being the number of all electrons in all atoms in unit volume, e and m 
are the electron charge and mass, respectively. Let C be a closed contour in 
the UJ upper half-plane defined as C = Ci U C2 , where Ci is the path following 
the real axis, except it goes around the origin w = in a small semicircle 
whose radius e tends to zero and C2 is a semicircle at infinity (see Fig. 1). If 
we denote by {ujj} all poles of the function /(w) in the region bounded by 
the contour C, the procedure described in Ref. [26] implies 



^Im/(.) = ./(0) + 2.^^H£(Zli^ 



(3.7) 



where Res means the residue. The dielectric function e{uj) has no zeros in 
the upper half-plane, so the only source of the poles of f{uj) inside C is the 
function g{uj) given by (|3.2p . Since g{uj) can be expanded in uj as follows |31j 



(3huj 



= 1 



2uj' 



UJ^ 



(3.8) 
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it exhibits in the upper half-plane an infinite sequence of simple poles at the 
(imaginary) Matsubara frequencies 

ujj = Res(g, ujj) = lVj [j = 1,2,...). (3.9) 

We thus conclude that the second moment of the static charge-charge density 
correlation function (|3.5p is expressible as 



dVr2/3(p(r)p(0))* 



1 

2^ 



1 

W) 



(3.10) 



In the u upper half-plane, the dielectric function e{Lo) takes real values only 
just on the imaginary axis and this fact ensures that the second moment 
p.lOp is real. Moreover, e{uj) decreases on the imaginary axis monotonically 
from eo > 1 (for dielectrics) or from oo (for conductors) at a; = iO to 1 at 
ijj = ioo; the convergence of the sum on the r.h.s. of (|3.10p is ensured by the 
asymptotic behavior (|3.6p . 

The representation (|3.10p is very useful in the high-temperature (classical) 
limit ptiLOp when each of the frequencies is much larger than 

LUp, the corresponding terms in the sum on the r.h.s. of Eq. (|3.10p vanish and 
therefore 



2,71 



1 

W) 



- 1 



(3.11) 



In the case of conductors, setting e(0) ioo implies the universal Stilhnger- 
Lovett sum rule [5]. In the case of a dielectric of static dielectric constant 
e(0) = eo, the above nonuniversal result has already been obtained by Chan- 
dler [33]; see also the article [34]. It should be noted that, at finite temper- 
ature, a quantum system always contains a finite fraction of ionized charges 
and thus, in a strict sense, it is conductor. In the case of the so-called di- 
electrics, the corresponding screening length is huge, much larger than the 
size of considered macroscopic samples, and therefore one does observe di- 
electric properties at experimental length scales. This point is well illustrated 
by exact results for atomic hydrogen I35j. 

The formula (|3.1ip can be rederived for dielectrics in an alternative way 
by extending screening and linear response ideas in classical Coulomb fluids 
[32ll6] to dielectrics as follows [36]. If we put an infinitesimal charge 5e at 
position r, the dielectric - test charge interaction Hamiltonian reads i?int = 
5e0(r), where 0(r) is the potential at point r created purely by the dielectric. 
By linear response theory, the average of this potential at some point r' is 
changed by the presence of 5e by 



W))se = -/3(i?int<^(r')) = ~/35e(0(r)0(r')), (3.12) 

where we suppose that ((/>(r)) = 0. The basic physical assumption about 
screening in dielectric media is that only a fraction (1 — I/eq) of the external 
charge is screened. Consequently, 5e surrounds itself with a polarization cloud 
of microscopic size carrying charge —(1 — l/tQ)5e. For |r — r'| much larger 
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than any microscopic scale, the l.h.s. of ()3.12p is thus the potential at r' due 
to a pointhke charge (1/eo — l)Se at r, i.e. 

The charge density p{r) is related to the potential (j){r) via the Poisson equa- 
tion A(p — —4:TTp. Taking first the Laplacian on r and then on r' in both sides 
of Eq. p.l3p . we obtain 

P{p{r)p{r')) = ^ (^1 - 1^ AS{r r'). (3.14) 

The second-moment sum rule (|3.1ip follows from this expression after an 
integration by parts, 

J d(r - r') (r ~ r' f (3 {p{v) p{y' )) = A (^1 _ j ^^"^ ~ ""'^ ^^"^ ~ ""'^ 

^ ^ ' - iV (3.15) 



27r \eo 

The formula (j3.10p provides a split of the integral value onto its classical 
and purely quantum-mechanical parts. In the quantum regime, it can be used 
to generate a systematic semiclassical expansion of the second charge-charge 
moment in powers of j3fi. In general, the r.h.s. of p.lOp is not expressible 
in terms of elementary functions. Perhaps the only exception is the jellium 
model of conductors whose dielectric function is adequately described, in 
the long- wavelength limit q — > 0, by the Drude formula with the dissipation 
constant taken as positive infinitesimal [29], 

£(^)-l- , I . V (3.16) 

LjJ{u} + 11]) 

Inserting (|3.16p with 77 = into (|3.10p and using the analog of the summation 
formula (fXS)l 

y ^ ^(acotha-l) (3.17) 

for a = we obtain 

i J A\r^m^)mY = -^5K)- (3-18) 

This result coincides with the microscopic finding [7,8J. Of course, the result 
p.lSp could be derived much simpler by inserting the Drude formula for the 
jellium (|3.16p directly into (j3.5p and then calculating the ensuing integral by 
using the Wcicrstrass theorem, see Eq. p. 191) which follows. The importance 
of the formula (|3.10p consists in the fact that it provides the division of the 
integral onto its classical and purely quantum-mechanical parts for all media. 
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3.2 Time-dependent case 

In the case of the jelHum model with the dielectric function (|3.16p . the Weier- 
strass theorem 

lim -^-r- ^ ( - I T i7r(5(a;) (3.19) 
,;-,o+ x±n] \x J 

{V denotes the Cauchy principal value) implies 

Im^ = -^sgn(c^) [Siu; - Up) + 5{uj + ujp)] . (3.20) 
Substituting this expression into Eq. (j3.4p . we arrive at 

i j A\r^P{p{t, r)p(0, 0))^ = "^^K) cos(wpi)- (3-21) 

The time-dependent second-moment sum rule for the unsymmetrized 
charge-charge density correlation function has been obtained microscopically 
in Refs. PITU]. 

i /dVr2/3(p(t,r)p(0,O))= ^ ^^'"^ 



27r 2 \ el'^^'^^p - 1 1 - e-P^^^p 

(3.22) 

It is easy to verify that the symmetrization of this expression leads directly 
to the result ([3?2T|) . 

The general formula (|3.4p is applicable to any conducting or dielectric 
medium, but since no other than jellium microscopic results are available we 
shall omit its analysis. 



4 Charge- current density correlations 

As concerns the charge-current density correlation function, using the spec- 
tral representations (|2.14p and (|2.18p we obtain 



(3oj 
(4^ 

9i^) 



2tt 



Im 



(AjAk) 



(4.1) 



As a check of the validity of this relation, we first multiply its both sides 
by Qk , then sum over k and finally use the spectral version of the continuity 
equation (|2.19p to arrive at (|3.1[) . 

The inverse Fourier transform of (|4.ip leads to two conditions: 







(4.2) 



12 



and 

/ d\xi(3{p{t,r)M0,O)y = Sh I dio exp(-iL^t)^Ini J-^. (4.3) 

Note that the sum rule (|4.3p fixes the correlation function of the fcth current 
component at time and the corresponding /th component of the total dipole 
moment at time t. The conditions (|4.2p and (|4.3p have the same form in both 
nonretarded and retarded regimes. 

The static t = Q version of the sum rule (|4.3p is trivial. Since g{Lo)lv[\ e^^{Ld) 
is an odd function of tj, we have 

J d\xip{p{r)MO)r ^0. (4.4) 

Let us now consider the time-dependent version of the condition (|4.3p for 
the jellium. Substituting the formula (|3.20p into (|4.3p . we get 

J d'rxip{p{t,r)M0,O)y = Ski^^^^smiio.t). (4.5) 

The same result has been obtained, after the symmctrization, by microscopic 
approaches [TUIITT] . 



5 Current- current density correlations 

For the current-current density correlation function, using the spectral rep- 
resentation (|2.18p we obtain after some algebra that 



/3(jfeji)t^.q = 7, 

ZTT 



fe-^)lm|ii^M^^j + ^Im^ 



(5.1) 

As a check, multiplying both sides of this relation by qi , summing over / and 
using the spectral version of the continuity equation (|2.19l) we arrive at (|4.ip . 

The Fourier component (|5.ip has a well defined value at q = since the 
two "dangerous" terms of type quqi/q^ just cancel with one another at this 
point, in both quantum and classical regimes. As a result, the current-current 
density correlation function is integrable in space and its zeroth moment reads 

di^ eM-'^ojt)^^lm— (5.2) 
-oo 47r^ e(^) 

The spatial integrability of the current-current density correlation function 
observed here is, specifically in the quantum regime, in contradiction with 
findings of microscopic theories for particles interacting through Coulomb 
forces only [7l[9lll0j: a discussion about this subject will be given below in 
Sect. 6. 
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5.1 Static case 



In the static t — case, the formula (15.21) can be rewritten in the form 



dV/30fe(r)jKO))^ = -4i 



duj 



Im/H, 



(5.3) 



where 



47r2 



6(C.) 



(5.4) 



the subtraction of the real sum + Wp from uP' jeito) has no effect due to the 
operation Im. The /-function possesses the symmetry /*(w) = /(— w) and 
vanishes at cx) due to the asymptotic relation (13. 6|) . We can therefore 

apply the integration techniques together with the general analytic properties 
of dielectric functions in the complex frequency upper half-plane presented 
in the second paragraph of Sect. 3.1, see Eqs. ()3.7|) - ()3.10|) . The final result 
reads 



/ 



d'r/3(ife(r)jz(0))^ 



2 - oo 

47r 27r ^ 



(5.5) 



where = 2Tij /{(3h) are the (real) Matsubara frequencies. 

In the high-temperature (classical) limit phtOp —^ 0, each of the frequen- 
cies is much larger than ujp, the corresponding terms in the sum on 
the r.h.s. of Eq. (j5.5p vanish and thus, for all conductor and dielectric media, 
it holds 



dV/3(jfe(r)j,(0))ei = fe^. 



(5.6) 



This classical result can be obtained in a simpler alternative way. In the 
classical regime, the position-dependent vector potential of the EM field can 
be eliminated from the momentum part of the Gibbs measure as well as from 
the definition of the current density by the Bohr-van Leeuwen theorem |37[ 
155] . Consequently, 



/3(jfe(r)jKr'))ci-/3e' 



N 



i=i j=i I 



(5.7) 



The momentum part of the one-particle Boltzmann factor is proportional to 
exp(— /3mt|2/2), so that (w^)ci = l/(/3"i). With regard to the definition of the 
plasma frequency p.6p we conclude that 



(5.8) 



This result has been obtained for the classical jellium longtime ago by Felder- 
hof ^i3j. The classical sum rule (|5.6p is a trivial consequence of the exact 
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relation ()5.8p . Moreover, higher moments of the current-current density cor- 
relation function vanish in the classical limit, 

y"dVr2^/30fc(r)ji(O))ei = 0, i=l,2,.... (5.9) 

In the quantum regime, the formula (|5.5[) can be used to generate a sys- 
tematic semiclassical expansion of the zeroth moment of the current-current 
density correlation in powers of (3h. For the jellium conductor with the Drude 
dielectric function (I3.16p . using the summation formula (I3.17P we obtain ex- 
plicitly 

J d\ P{jk{r)ji{0)r = 4i^g(c^p). (5.10) 

This result is not confirmed by microscopic theories for particles interacting 
only through Coulomb forces [TllQllTO] which predict the spatial nonintegrabil- 
ity of the current-current density correlation function in the quantum regime 
(see Sect. 6). 



5.2 Time-dependent case 

In the special case of the jellium, using the relation (|3.20p in Eq. (|5.2p . the 
generalization of the static formula (|5.10p to time difference t is found to be 



dV/3(jfe(t,r)j,(0,O))^ = Ski-^giiOp)cos{copt). (5.11) 

47r 



6 Problematic points in current-current fluctuations 

There are certain problematic points in the results for the current-current 
density correlation function, in the fluctuational theory presented here as 
well as in microscopic approaches [711^1 lOj. 

We start with an obvious inconsistency having the origin in the formula 
(|5.ip . We shall consider the jellium with the dielectric constant (|3.16p . Using 
the relation p.20p and the similar one 



with a modified plasma frequency ujp — ^ujp + {cqY, the partial inverse 
Fourier transformation of (|5.ip in time leads to 



w2 



—giujp) cos{ujpt) \du — 1 + g{ujp)cos(u}pt)^^ 



(6.2) 



15 



In the static t = case and in the classical limit, which corresponds to setting 
giiOp) = g{ijjp) = 1, this formula is equivalent to 



d3re"''i-/3(jfe(r)jKO))c 



p 
An 



( r qkqi \ , Qkqi 



{cqf 



2 / ' „2 



= '^'^'^ + -^il^qi - hiq^) + 0{q^), (6.3) 

which is an infinite series in the wave- vector components. On the other hand, 
as follows from the exact analysis between Eqs. (|5.6p - (|5.9p . only the ze- 
roth moment of the classical current-current density correlation function is 
nonzero, so that it must hold 

d\e-'^-'^f3{Mr)ji{0)U = 4;^. (6.4) 

Our expression ()6.3p coincides with this exact formula only in the leading 
absolute term which was in the center of interest in Sect. 5. The erroneous 
presence of higher-order terms is intuitively related to the fact that Rytov's 
fluctuational theory works with the local dielectric function e{uj) which is the 
long-wavelength limit of the q-dependent one e(uj, q), e(a;) = limq^o e(w, q). 
A complete theory should involve the whole function e{oj, q). We can antici- 
pate the validity of Rytov's fluctuational theory only in the leading q order, 
in which the q-dependence of e(a;, q) can be ignored without any approxima- 
tion. In fact, without saying it, this strategy has already been adopted. 

In microscopic approaches for purely Coulomb systems [7jl9jfT0]. the ob- 
tained result for the (partial) Fourier transform of the current-current density 
correlation function 



qkqi\ , / ^ / ..Qkqi 

—p- I + g(LJp)cos(ujpt)—p- 



(6.5) 



differs from the present one (|6.2p in the prefactor to the transverse component 
which does not depend neither on time t nor on f3h. The expression (16. 5p 
works well in the static {t — 0) classical {/STi O,g{ojp) = 1) regime, where 
it reproduces the exact formula (|6.4|) . As soon as the time difference t ^ 
or the regime is quantum {f3h ^ 0,g{ujp) > 1), the two "dangerous" terms of 
type quqi/q^ do not cancel with one another in (j6.5p . Since 

d'g jc^-rikqi _ 9' [ d^g ^^^.^ 1^ _ 5^ 1 ^ ^^ ^^ 



{2'itY q^ dxkdxi J (27r)3 q"^ dxkdxi Airr' 

the current-current density correlation function has a nonintegrable algebraic 
decay at large distances. 

The discrepancy between our (retarded) and microscopic (non-retarded) 
results for the current-current density correlations is probably analogous to 
that for the induced electric potential and field correlation functions [T5ll2T] . 
Our analysis does not exclude an algebraic long-range behavior of the current- 
current density correlation function, but this decay must be fast enough to 
ensure the spatial integrability of the correlation function. 
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7 Classical limit of current-current fluctuations 

The classical limit of the current-current density correlation function for the 
jellium was investigated longtime ago by Felderhof, in the nonrctarded regime 
[2] as well as for the retarded plasma fully coupled to the classical radiation 
[3]. Changing the notation k ^ q and using the substitution s = mu/q, his 
final formula (13.5) in Ref. [3] can be rewritten after some straightforward 
algebra as follows 



where 



'2^ 



q'^ J \ Ziq'^/{muj) 



Z^{uj,q) = l-ujp dx 



{uj + vq/2 — x)'- 



(7.1) 
(7.2) 



and 



Zi{w,q) = ( — 



17; 



-6aix) 



(7.3) 



with 77 — > 0+; the function 6a{x) and the parameter a are defined by 

^,(x)=.^e->^ a^q.[^. (7.4) 
V7ra V P"^ 

In the limit q ^ (a ^ 0), 5a{x) becomes one of the standard definitions 
of the (5-function. Consequently, 

Zo{uj,q) ^ e{iu), Z,{iu,q)(^Y ^ e{u;)~(^)\ (7.5) 

9^0 KmUoJ q^O \ UJ / 

where e(w) is the dielectric function of the jellium (|3.16p . Inserting these ex- 
pressions into (|7.ip and comparing the resulting formula with the classical 
version {g{uj) = 1) of the current-current density correlation function (15. ip . 
we see the equivalence of the two approaches in the considered limit q — > 0. 
However, when q is finite, the approaches give different results. Felderhof's 
classical formula turns out to be correct for an arbitrary value of q since it 
leads to the exact relation (|5.8p or its equivalent (|6.4p 3 , . From the discussion 
in Sect. 6 we know that this is no longer true for the classical version of our 
formula (j5.ip . This fact might be a motivation for searching new improve- 
ments of Rytov's fluctuational theory which reproduce classical results for 
any value of q and simultaneously describe adequately correlation functions 
in the quantum regime. 
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8 Conclusion 

In this paper, by using Rytov's fluctuational theory we have derived a se- 
ries of known or new sum rules for the bulk correlation functions of charge 
and current densities in conductor and dielectric media, fully coupled to the 
radiation field (retarded regime). Few important facts are worth to mention. 

In the Fourier q space, we can anticipate the validity of Rytov's macro- 
scopic approach only in the leading q order. In this order, the q-dependence 
of the dielectric function e{Lu,q) can be ignored without any approximation. 

In the case of the static t — sum rules, an analysis in the complex 
frequency upper half-plane was used to split these sum rules onto their clas- 
sical and purely quantum-mechanical parts. The classical part coincides, for 
both conductors and dielectrics, with the results obtained previously by other 
methods. The quantum-mechanical part was checked on the jellium model of 
conductor. 

As concerns the time-dependent sum rules, their validity was tested on the 
jellium, too. The retarded results for the charge-charge and charge-current 
density correlation functions coincide with the non-retarded ones obtained 
previously by microscopic approaches to the jellium model. The general for- 
mulas are applicable to any media, but no other than jellium microscopic 
results are available to make a comparison. 

The most interesting quantity is the current-current density correlation 
function. Within microscopic theories in the non-retarded regime, this cor- 
relation function turns out to be spatially integrable in the classical limit, 
but it becomes nonintegrable in the quantum case due to an algebraic decay 
at large distances. This contradicts our retarded prediction (|5.2|) about 
a finite value of the integral over the space, in both classical and quantum 
cases. 

A generalization of the present results to inhomogeneous situations, like 
a surface contact between two distinct media, is left for the future. 
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